Mapping broadband single-photon wavepackets into an atomic memory 
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We analyze a quantum optical memory based on the off-resonant Raman interaction of a single 
broadband photon, copropagating with a classical control pulse, with an atomic ensemble. The 
conditions under which the memory can perform optimally are found, by means of a 'universal' mode 
decomposition. This enables the memory efficiency to be specified in terms of a single parameter, 
and the control field pulse shape to be determined via a simple nonlinear scaling. We apply the 
same decomposition to determine the optimal configurations for read-out. 
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Introduction The conversion of optical information 
between photons and atoms forms the basis of a quantum 
interface that is a critical component of quantum commu- 
nications networks and distributed quantum computers 
U 0, 9 • Such quantum information processing schemes 
require the ability both to move quantum information 
between nodes of a network, and to store the informa- 
tion. An important class of quantum optical memories 
is based on the interaction of individual photons with 
atomic ensembles, in which the information is transferred 
coherently from a single photon to a collective excitation 
of the atoms. If active feedback Q is not used, these 
memories are generally optically dense, with the density 
controlled dynamically using an ancillary field [H, H, @] . 

In this paper we analyze a prototypical memory based 
on the off-resonant Raman interaction of a classical 
pulsed control field and a broadband signal photon in 
an atomic medium (see Fig. [!}. The temporal struc- 
ture of the signal photon is transferred by the control 
to a long-lived collective atomic excitation [5(, or spin 
wave. This differs qualitatively from previous narrow- 
band schemes 8] which demonstrated that quadrature 
squeezing could be transferred from optical fields to a col- 
lective atomic spin. We show how proper shaping of the 
control field allows the mapping of an input wavepacket 
of arbitrary temporal shape to an output wavepacket of a 
potentially different temporal shape. The dynamics are 
closely related to those investigated in proposals for en- 
tanglement generation via spontaneous @, ul and, more 
recently, stimulated 1(| 11] Stokes scattering. However 
the photon storage process is distinct from these, in that 
it exhibits an explicit time reversal symmetry; evinced by 
its fundamental mode structure. Spontaneous emission 
is suppressed as long as the excited state remains empty. 
Lossless unitary storage of broadband single photons - 
such as are commonly used in cryptographic and telepor- 
tation experiments [H, [HJ _ can therefore be implemented 
by detuning sufficiently from resonance 17). Departure 
from resonance makes our scheme robust against inho- 
mogeneities in the ensemble, so that solid state absorbers 
(e.g. semiconductor charge quantum dots [HI]) could be 
substituted for the atoms. In addition, changing the de- 



tuning of the control pulse between storage and retrieval 
allows for control over the frequency of the output state. 
In the following we consider propagation in the one di- 
mensional limit; a fully three dimensional model will be 
considered elsewhere. 
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FIG. 1: Left: the level structure of the i atom comprising 
a quantum memory for broadband photons, with bandwidth 
8. Right: a schematic of the read-in process for the quantum 
memory. 



Model The signal and control fields are Raman res- 
onant, with center frequencies u> s , lu c , respectively. The 
classical control at time t and position z is represented 
by the Rabi frequency f2(r), where r = t — z/c is the local 
time. The signal and spin wave amplitudes are described 
by the slowly varying annihilation operators A(t, z), and 
B(t,z), respectively (lij |. The spin wave is a collective 
coherence of the form B(t,z) cx Y,p |l/3)(3 /3 |e- i ^ c " ; <=) r , 
where the index [3 runs over all atoms with position z 
[Toj ] . If the common detuning A of the signal and control 
pulses from single photon resonance is much larger than 
the signal bandwidth S, the control Rabi frequency Q, 
and the control bandwidth, the excited state |m,-) can be 
adiabatically eliminated. If the ensemble is prepared in 
the collective state |0) = |L), and if the population 
of the metastable state |3j) is assumed to remain negli- 
gible, a linear theory can be used. The Maxwell-Bloch 
equations, in the slowly varying envelope approximation, 
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are then found to be 

[d T -i\n(T)\ 2 /r]B(T,z) = 
[d z -i\ K \ 2 /r}A(r,z) = 



- K *n*( T )A(T,z)/T; (1) 
K n(r)B(T,z)/T, (2) 



2 



where k is the signal field coupling and Y = A — 17 = 
|r|e _l9 is the complex detuning, with real phase 9. 7 
arises from dephasing processes, including spontaneous 
emission. We have not included the Langevin noise oper- 
ator which formally accompanies these loss terms, since 
its contribution vanishes when normally ordered expec- 
tation values of A and B are taken. We neglect the slow 
decay of the spin wave over the short timescale of the 
memory interaction. 

We now introduce a new set of scaled coordinates: 
the memory time e(r) = Cuj(t)/uj(T), where T is 
the duration of the interaction, and the effective dis- 
tance C,(z) = Cz/L, with L the length of the ensem- 
ble. Here lo{t) = JJ" |f2(r')| 2 dr' is the integrated Rabi 
frequency and C = \k\ Llo(T)/\Y\ is a coupling pa- 
rameter. We define dimcnsionlcss annihilation operators 
a(e,C) = ^Lo{T)/Ce~ ix ^A(T, z)/&{t) for the optical 
field (assuming k is real for simplicity) and 0(e, £) = 
^L/Ce~ lx ( T ' z ^B(T, z) for the spin wave, where the ex- 
ponent is x(t, z) = [u>(t) + |K| 2 z]/r. The first term in \ 
describes a Stark shift due to the control field; the sec- 
ond represents a modification of the signal group veloc- 
ity. With these changes, the equations of motion reduce 
to the simple coupled system d^a = e l9 /3; d e (3 = —e l9 a. 
The solution of these equations then holds for all control 
pulse shapes and arbitrary inputs. The coupling parame- 
ter C sets the size of the region in (e, C)-space over which 
the memory interaction is driven. For the case consid- 
ered here, an atomic ensemble, C can be re-written in 
the form C = (irafh/me) 1 / 2 f^/NJ^/ '(\Y\A), where / 
is the geometric mean of the oscillator strengths for the 
signal and control transitions, A is the cross-sectional 
area of the control field, and N a (N c ) is the number of 
atoms (photons) interacting with (comprising) the con- 
trol pulse. Here a/ is the fine structure constant, and m e 
is the electron mass. 

The memory read-in and read-out must be unitary to 
function correctly. We now show that canonical evolu- 
tion of the field operators a and /3 is guaranteed by the 
classical structure of Eqs. (|l|2p in the dispersive limit 
A ^ 7. In this case, the phase 9 vanishes and the follow- 
ing continuity relation is satisfied: d^a^a + d c j3^ (3 = 0. 
Integration of this expression over a square in (e, £)-spacc 
yields the flux-excitation conservation condition 



N a (C) + Nf)(C) = N a (0) + N (O), 



(3) 



where the number 



operators 

/ C at( ejC ) a ( e;C ) de , and Np(e) = &(e, C)/?(e, C) dC, 
count the number of signal photons at an effective 
distance £, and the number of excitations of the spin 
wave at memory time e, respectively. Eq. ([3]) must 
hold for arbitrary initial amplitudes 00(e) = a(e,0) 
and /3o(C) = /3(0>C)> which fixes the transformation 
{00(c), ^ (C)} -> {«c(e),/3c(C)} as unitary (where 
oc(e) = o(e, C) is the signal amplitude at the exit face 



of the ensemble, and /3c (C) = ls the spin wave 

amplitude at the end of the read- in process). This allows 
the dynamics to be decomposed into a set of independent 
transformations between light-field and spin wave modes 
TH 14 1 . In what follows we therefore concentrate on the 
dispersive limit, and consider the case r — » A; 9 — * 0. 
The solution of the dynamical equations is expressed by 
the scattering relations 13, 15] 



a c (e) = f [die - x, C)a {x) + G (C - x, e)/3 (x)]dx; (4) 
Jo 

/?c(C)= / [G 1 (C-x,C)f3 (x)-Go(C-x ) Oao(x)]dx.(5) 
Jo 

The integral kernels are given by Go(p,q) = Jo (2y/pq), 
and Gi(p, q) = 6{p) - 0(p) Ji (2y/pq) y/q/p, with the n th 
Bessel function of the first kind denoted by J„, and where 
the Heaviside step function ensures that the convolu- 
tions in Eqs. (|4I5[) respect causality. 

The integral kernels Go.i - as they appear in Eqs. (1415 j) 
- share symmetry under reflection about the line C — x = 
y, where y stands for the independent variable: either 
e or £. This symmetry, along with Eq. ([3]), allows us 
to decompose the kernels using input and output modes 
related by time reversal (or equivalently space reversal) 
as follows: 



G (C 
Gi(C- 



(6) 
(7) 



where {4>i} is a complete orthonormal set of real mode- 
functions and where the real, positive singular values sat- 
isfy the constraint Xf + fif — 1, Vi. 

Read-in The ensemble begins the read-in process in 
the state |0) , and we are free to replace /?o(C) by its expec- 
tation value; /?o(C) ~^ 0. With the above decomposition, 
Eq. ([5]) then describes a mapping of the optical input 
mode <pi(C — e) to the output spin wave mode 4>i(C), with 
transfer amplitude — Aj, for each i. Transforming from 
memory time e back to local time t, the normalized input 
modes are written $;(r) = ^C/w(T)e ix ( T '°)f2(T)</>j [C - 
e(r)]. The read-in efficiency can be quantified by evaluat- 
ing the expectation value of the spin wave number opera- 
tor (Np(C)) at the end of the read- in process. Expanding 
an incident signal wavepacket £(r) using the the read- 
in efficiency is expressed as (Np(C)) = J2i ^i\£.i\ 2 > with 
the i th overlap given by & = J Q ^*(t)$j(t) dr. When 
(N@(C)) — 1, the read-in works perfectly. In Fig. [2] the 
first five transfer amplitudes are plotted as a function of 
C. These are found using the eigenvalue equation 



Jo{2y/xy)^i(y) dy = A^x), 



(8) 



which we solve numerically using a 500 by 500 square 
grid. It is desirable to limit the energy of the control 



FIG. 2: The five largest singular values of the kernel Go, 
plotted as a function of the coupling parameter C. 




FIG. 3: Left: the intensity {A t (t, z)A(r, z)} of a Gaus- 
sian signal photon with wavepacket amplitude £(r) cx 
exp {—2 In 2[(t — to)/<t] 2 }, as it propagates through an atomic 
ensemble with C = 2. Here a = T/8, r = 2T/3. Right: the 
optimized control field intensity is shown alongside the initial 
signal field intensity (scaled for clarity). 



pulse, so we should find the minimum coupling parame- 
ter C which permits complete storage of the signal. For 
C « 2 the lowest mode achieves its optimal efficiency 
Ai ~ 1, but higher modes remain poorly coupled. The ef- 
ficiency of the memory is therefore maximized by setting 
£ = l; = 0, so that (N P (C)) = A? w 1 (for C > 2). 
To do this, it is necessary to shape the control field so 
that $i(t) = £(t). If Gaussian optics (.4 ~ cL/uj s ) are 
used to illuminate a region a few cm long in a typical 
atomic vapour (/ ~ 1) of modest density (~ 10 20 m~ 3 ), 
with 100 nJ control pulses, a 1 ps photon wavepacket 
can be stored optimally, with C — 2. In practice mode- 
matching could be achieved through measurement and 
feedback: the signal and control field sources are locked 
to a phase reference and operated in pulsed mode. The 
control pulse profile is characterized [la ] , and augmented 
until the transmission of the signal is minimized. Figure 
[3] shows the result of a simple optimization to find the 
control pulse shape which modematches the lowest input 
mode to a Gaussian signal photon, for C = 2. The pho- 
ton is absorbed with (Ng(C)) « 0.96; the small transmis- 
sion probability is due to the limitations of our numerical 
modematching optimization. 

Read-out Once a properly modematched photon has 
been read in to the quantum memory, the ensemble is 





FIG. 4: The photon retrieval probability Af, for forward read- 
out, plotted as a function of the read-in and read-out coupling 
parameters C and C r . 



left in the output mode </>i(C), with probability ampli- 
tude — Ai . We now consider the effect of sending a second 
control pulse, propagating in the same direction as the 
initial control pulse, into the ensemble. The center fre- 
quency, bandwidth, and intensity of this read-out pulse 
may differ from that of the first control pulse (herein the 
read-in pulse). Let us use a superscript r to indicate 
those quantities associated with the read-out. We ne- 
glect decoherence and dephasing of the spin wave over 
the storage period and set B r (0, z) = B(T, z). This pro- 
vides us with one boundary condition; the second is that 
the signal field begins in its vacuum state at the start of 
the read-out process, (</V£(0)) = 0. The efficiency of the 
read-out depends upon the degree to which the spin wave 
mode ipi(z) = y/ C / Le lx ^ T,z ^ (f>\{C z / L) (written in terms 
of the ordinary spatial variable z) overlaps with the in- 
put modes = yJC r /Le^'^tf [C r (l - z/L)\ for 
the read-out process. The functions {</>[} solve the eigen- 
value equation ©, with C replaced by C r . A measure of 
the efficiency of the memory is the expectation value of 
the output photon number operator Af = (N^(C r )) = 
Si K 2 \fi\ 2 ' with the read-out overlaps defined by 
fi = J ip* (z)^lj (z) dz. The parameter Af is the probabil- 
ity of retrieving a photon from the ensemble at read-out, 
given that a single modematched photon was sent in with 
the read-in pulse. If the detuning does not change too 
much, so that |A r - A|/A < A7(|k| 2 L) - y/N r c /N a , 
then the phases x, \ r °f the spin wave modes approxi- 
mately cancel, and then the stored spin wave is phase- 
matched to the readout modes. In Fig. [4] the variation 
of Af, under this approximation, is plotted as a function 
of the read-in and read-out coupling parameters C and 
C r . If C — C r , then the lowest read-out mode is just the 
mirror image of the spin wave mode. For small C, the 
spin wave mode is monotonic, and relatively flat; f± is 
therefore large. However the transfer amplitudes Ai and 
X[ remain small, so the retrieval probability is low. In- 
creasing Ai requires a larger C, but this produces a more 
asymmetric spin wave mode, and f\ falls. It is then nec- 
essary to increase C r above C, so that higher modes, with 
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which the spin wave mode overlaps significantly, are effi- 
ciently coupled to the optical field. The retrieval proba- 
bility is maximized along the line C«2, which represents 
the optimal coupling for the read-in process. However, a 
read-out coupling parameter in excess of 10 is required to 
achieve AT > 0.95. Note also that modematching to the 
lowest mode at read-in is the best strategy for maximiz- 
ing the memory efficiency. Modematching to a higher 
mode, or some combination of modes, simply increases 
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the optimal read-in coupling above C = 2 

The time reversal symmetry between the input and 
output modes makes the read-out for this scheme a non- 
trivial problem. Simply repeating the read-in process (so 
that C = C r ) results in poor performance of the mem- 
ory. However, the dramatic increase in coupling strength 
required to extract the stored excitation fully, may make 
a naive increase in control pulse energy at read-out pro- 
hibitively difficult to realize. An alternative method to 
boost the coupling is to reduce the bandwidth of the read- 
out pulse, along with its detuning A r . The photon recov- 
ered from the memory in this way would be frequency- 
shifted (according to the Raman resonance condition), 
and temporally stretched (since its bandwidth would be 
diminished as well). Such a memory would act as a 'pho- 
ton transducer', storing broadband photons and convert- 
ing them to narrowband photons with tunable frequency 
on demand. 

Note that switching the propagation direction of the 
read-out pulse sends cj>i(Cz/L) — > <f>i [C(l — z/L)], so 
that the spin wave mode overlaps exactly with the low- 
est read-out mode with C r = C, and we should obtain 
Af = Af . Unfortunately the read-out process is no longer 
phasematched in this situation, and the overlap integrals 
fi vanish. However, a solid-state implementation might 
allow this kind of reverse read-out with the use of quasi- 
phasematching, in which the sign of the read-out coupling 
parameter C r is periodically flipped along the length of 
the ensemble. 

The phasematching problem is obviated in the limit of 
vanishing Stokes shift, but then the ground state |0) must 
be prepared artificially with high purity: if the state |3j) 
is initially populated, or if selection rules allow residual 
coupling of the control to the ground state, the memory 
fidelity at the level of single quanta is greatly reduced. 
Furthermore, if signal and control are spectrally indis- 
tinguishable, another degree of freedom should be used 
to differentiate between them. Typically, polarization af- 
fords discrimination to a part in ~ 10 6 , but then the 
control should not contain more than a million photons. 
A large Stokes shift is therefore desirable, and correct 
phasematching at read-out is crucial for efficient retrieval. 
The above considerations demonstrate the importance of 
propagation effects in the design of an optical quantum 
memory. 

Summary We have shown that the off-resonant Ra- 
man configuration for a A-type ensemble quantum mem- 



ory can be used to implement deterministic, controllable, 
and unitary transfer of the temporal structure of broad- 
band single photons to a stationary spin wave, in the 
adiabatic regime. The dynamics are understood and op- 
timized using a universal mode decomposition, valid for 
all control pulses and arbitrary input states. The modes 
are computationally simple to evaluate, and only a few of 
them are required to approximate the interaction faith- 
fully. The optimal fidelity of the memory depends only 
upon a single dimcnsionless parameter (C), which defines 
an equivalence class for memories with different physical 
implementations. 
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